FERMIL AB-PUB-00 /035-T 
BNL-NT-00/14 
February 1, 2008 



Soft Double— Diffractive Higgs Production 
at Hadron Colliders 

Dmitri Kharzeev* a ) and Eugene Levin t b )' c ) 

°) Physics Department and RIKEN-BNL Research Center, 
Brookhaven National Laboratory, 
Upton, NY 11973 - 5000, USA 

b > HEP Department, School of Physics, 
Raymond and Beverly Sackler Faculty of Exact Science, 
Tel Aviv University, Tel Aviv 69978, ISRAEL 

c ) Theoretical Physics Department, 
Fermi National Accelerator Laboratory, 
P.O. Box 500, Batavia, Illinois 60510, USA 



Abstract 

We evaluate the non-perturbative contribution to the double-diffractive production of 
the Higgs boson, which arises due to the QCD scale anomaly if the mass of the Higgs Mjj 
is smaller than the mass of the top quark Mt, Mh < Mt- The cross section appears to 
be larger than expected from perturbative calculations; we find an = 0.019 -j- 0.14pbarn 
at the Tevatron energy, and an = 0.01 0.27pbarn at the energy of the LHC. 
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Figure 1: Lego-plot for double Pomeron Higgs production process. 

1 Introduction 

In this letter we suggest a new mechanism for "soft" double-diffractive production of Higgs 
boson. We consider three reactions 

p + p — > p+ [LRG] +H+ [LRG] +p ; (1) 
p + p — > Xi + [Li?G] + H + [LRG] + X 2 ; (2) 
p + p — > p + + if + [L_RG] + X 2 ; (3) 

where LRG denotes the large rapidity gap between produced particles and X corresponds to a 
system of hadrons with masses much smaller than the total energy. These reactions have such 
a clean signature for experimental search (see Fig.l, where the lego - plot is shown for reaction 
of Eq. (|])) that they have been the subject of continuing theoretical studies during this decade 
( see Refs.0, 0, §, §§,§). 

The main idea behind all calculations, starting from the Bialas-Landshoff paper [|TJ, is to 
describe the reactions of Eq. (|T|) and Eq. (0) as a double Pomeron (DP) Higgs production ( 
see Fig.2 ) . In Fig. 2, the Pomerons are the so-called "soft" Pomerons for which one uses 
the phenomenological Donnachie-Landshoff form ( see Ref. 0] ), while the vertex 7 can be 
calculated in perturbative QCD. 

We can demonstrate the problems and uncertainties of such kind of approach by considering 
the simplest pQCD diagram for the double Pomeron Higgs production ( DPHP ) (see Fig.3-a). 
This diagram leads to the amplitude 

M(qq - qHq) = 2 -2g H J d " Q± 4 a s (Q 2 ± ) (Qx,± ■ 4,±) , (4) 



2 




Figure 2: Double Pomeron Higgs production process. 

where gu is the Higgs coupling that has been evaluated in perturbative QCD ||. For the 
reaction of Eq. ([[]), \t\\ = \Q±_ — Qi,±| ~ \t%\ = \Q± — Q2.1I ~ 2/£? e ; and therefore, 

M (g + g-,g + iJ + g) oc / (5) 

Eq. (^|) has an infrared divergence which is regularized by the size of the colliding hadrons. In 
other words, one can see that already the simplest diagrams show that the DP Higgs production 
is, in a sense, a "soft" process. Taking into account the emission of extra gluons denoted in 
Fig.3-b as Pomeron builders, we recover the exchange of the "soft" Pomerons. 

Nevertheless, the emission vertex for the Higgs boson can still be calculated in pQCD since 
the typical distances inside the quark triangle in fig.3-a are rather short oc where My is 

the mass of t-quark. The coupling gjj has been evaluated in Ref . |3| and is given by 

g% = y/2G F a 2 s (M 2 H ) iV 2 /9vr 2 , (6) 

where N is a function of the ratio Mt/Mh which was calculated in Refs. || |3|. 

In this paper we consider an alternative approach to DPHP, in which we estimate the value 
of the cross section from non-perturbative QCD. In section 2 we review a non-perturbative 
method suggested by Shifman, Vainshtein and Zakharov || for the evaluation of the coupling 
of the Higgs boson to hadrons; it is valid if the mass of the Higgs is smaller than the mass of 
the top quark. In section 3 we develop a method of obtaining the DPHP cross section using the 
approach of Ref. |9| . The problem of survival of large rapidity gaps (LRG) will be discussed 
in section 4. We conclude in section 5 with discussion of our results and of the uncertainties 
inherent to our approach. 
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Figure 3: Double Pomeron Higgs production in the Born approximation (Fig. 3-a ) and in 
leading log approximation (Fig. 3-b ) of pQCD. 



2 The coupling of Higgs boson to hadrons in 
non-perturbative QCD 

To evaluate the non-perturbative coupling of the Higgs boson to hadrons, we need to have a 
closer look at the properties of the energy-momentum tensor in QCD. The trace of this tensor 
is given by 



a 

a — ~~2^~ a P 



Y m *(! + imMm + Y 

l=u,d,s h=c,b,t 



m h {l + j mh )QhQh, 



(7) 



where 7 m are the anomalous dimensions; in the following we will assume that the current quark 
masses are redefined as (1 + j m )m. The appearance of the scalar gluon operator in (|7|) is the 



consequence of scale anomaly JTOf, [11|. The QCD beta function can be written as 



P(g) = -b-^ + ..., & = 9-~7l fc , 
107T^ 3 



(8) 



where is the number of heavy flavors (c, b, ..). Since there is no valence heavy quarks inside 
light hadrons, at scales Q 2 < 4m| one expects decoupling of heavy flavors. This decoupling was 
consistently treated in the framework of the heavy-quark expansion |J; to order only 
the triangle graph with external gluon lines contributes. Explicit calculation shows B that the 
heavy-quark terms transform in the piece of the anomalous gluonic part of 0°: 



m hQhQh 

h 



2 9 



32vr2 G 



(9) 



It is immediate to see from (||), (|7]) and (H) that the heavy-quark terms indeed cancel the part 
of anomalous gluonic term associated with heavy flavors, so that the matrix element of the 
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energy-momentum tensor can be rewritten in the form 

Q a = M G ^a G a^ + £ mjgftj (1Q) 
" l=u,d,s 

where heavy quarks do not appear at all; the beta function in (3.10) includes the contributions 
of light flavors only: 

/%) = -9 I f^ + --- (11) 

Because the mass of the Higgs boson Mh is presumably large, its coupling to hadrons 
involves the knowledge of hadronic matrix elements at the scale Q 2 ~ M H , at which the heavy 
quarks in general are not expected to decouple. However, if the Higgs boson mass Mh is smaller 
than the mass of the top quark Mt, one can still perform expansion in the ratio Mh/Mt; we 
expect this to be a reasonable procedure if Mh < 100 GeV. In this case, one finds 

2 n 2 

M T tt -> y —G aPa G a nB + ... (12) 

Since the mass of the hadron is defined as the forward matrix element of the energy- 
momentum tensor, the expression Eq. (0) leads to the following Yukawa vertex for the coupling 
of a Higgs boson to the hadron: 

i i _ i i 2M 2 

2^Gl ■ H<f> 2 ■ (h\M T tt\h) = 2*Gl (13) 

this relation is valid in the chiral limit of massless light quarks (see Eq. (|TUD); Mt is the mass 
of the heavy quark and M is the hadron mass. We put the number of light quarks Np = 3 
and the number of colors N c = 3; <ph and H are hadron and Higgs operators. Note that, as a 
consequence of scale anomaly, Eq. (|13|) does not have an explicit dependence on the coupling 

3 Estimates for double Pomeron Higgs production cross 
sections 

3.1 General formulae for double Pomeron Higgs production 

The amplitude for Higgs production in the Pomeron approach is given by ( see for example 
Refsj,0) 

/ S \ a H*2) / s \ ap(*x) 

M(h+h -»• h+H+h) = g x {h) ■ g 2 {t 2 ) ■ i{t u t 2 ) ■ rj+^rj+ih) ■ —J . (_j ,(14) 

where S\ = (P\ + qi) 2 and s 2 = {P 2 + q\) 2 ( P\ t2 are momenta of incoming hadrons ); ?7+(t«) is 
a signature factor, which for the Pomeron is 

V+ (U) = i + tan-'l 7 ^^) , (15) 
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where ap(t) is the Pomeron trajectory, ap(t) = 1 + Ap + a' P t, with Ap ps 0.08 0; all other 
notations are evident from Fig.2. 

The cross section for DPHP in the central rapidity region (yjj = 0, where is the rapidity 
of the produced Higgs boson) can be written down as 



da 

dyHdhdh'""- = { ' 

where P[ are momenta of recoil hadrons, while pn is the momentum of the produced Higgs 
boson. 

Performing all integrations and recalling that s\ ■ S2 — M\ ■ s we obtain 



do 2g 2 (t 1 )-g 2 {t 2 )- 1 2 {t l M) ( _s_ 

dyndhd^"- ~ 7r(167r) 2 ' \M 2 H t 



2 A P 

a' ln(s/M%) [ti+t 2 ] 



e a' p ]n(s/M^)[ti+t2\ _ ^) 



We will assume that 7(^1,^2) is a smooth function of £1 and £2 in comparison with gi{t\) 
and #2 (£2)- Indeed, the t-dependence of is related to the quark distribution inside the hadron 
while the t-dependence of 7 is determined by the mean transverse of gluon inside the Pomeron. 
The typical scale for this momentum is l/ctp ~ 4 GeV 2 which is much larger than the typical 
momentum of a quark in a hadron ( sa 0.1 GeV 2 ). 

Using this assumption together with the simplest Gaussian parameterization for the vertex 
9i(ti) = <7i(0) exp(— Rq \U\) we obtain 

\VH=0 = Ua D /„/»/2M2 • 7 (*1 = = 0) 



cW" H ^[levr^^/M^)] 2 ,w " 7 V M lr 



Recalling now the well-known relation between the total and elastic cross sections for the one 
Pomeron exchange, namely, 

= fiiW = £«1 c a) a 7 (19) 



cr fot (s) 167r5 e /(s) V s 
where B d = 4R^ + 2a' P Ins, one can derive 

^\, B =o = 7 2 ( (l = 0, (2 = 0)x^x^(^ So ) . (20) 

There is only one unknown factor in Eq. (p0[) , namely, 7 2 (ti = 0, t 2 = 0). In the next subsection 
we present the estimates for this factor using the non-perturbative approach that has been 
discussed in the section 2. 
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Figure 4: Higgs emission from glueball. 

3.2 The production vertex j(ti = 0,^2 = 0) 
Our estimate of 7(ti = 0, £2 = 0) consists of two steps: 

1. For positive values of h = t 2 = m 2 glueball we can obtain 7^ = m 2 glueball , t 2 = m 2 glueball ) from 
Eq. (PD; 



2. Using Eq. (|I"4]) we can make the analytic continuation to the region ti < and £2 < 0, 
which corresponds to the scattering process. 

We will assume that there exists a tensor 2 ++ glueball which lies on the Pomeron trajectory, 
namely, that its mass satisfies the following relation: 

Mt = m 2 glueball ) = 1 + A + a' P (0)m 2 glueball = 2. (21) 

There is no undisputed experimental evidence for such a meson but lattice calculations give 
for its mass m g i ue baii = 2.4 GeV [T3[]. This mass is a little bit higher than can be expected 



from Eq. (|2T|) with the experimental a' P (0) = 0.25 GeV' 2 0. On the other hand, it is possible 
to describe experimental data using a smaller value of a' P (0) ~ 0.17 GeV~ 2 which is needed 
to satisfy Eq. (|2~ID with m g i ue b a ii = 2.4 GeV, assuming the presence of substantial shadowing 
corrections [jUR . 

For the diagram in Fig. 4 the vertex ■jgiuebaii can be easily evaluated from Eq. (0); it is equal 

to 

glueball = 2* Gl 2m |p" . (22) 
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One can see that Eq. (|H|) leads to the contribution described by Fig. 4. Indeed, for U 



2 

m glueball 



ri+(ti) — -r • (23) 

7r ct P ( rn glueball ■ ■ t, 



(A more detailed discussion of the analytic properties of the Reggeon exchange can be found 
in Ref. p5| ). Using Eq. (|23| ) and comparing Eq. ( |T4| ) with the diagram of Fig.4, we conclude 
that 



7T 

l(h = ml lueball ,t 2 = m 2 gluebaU ) = - a' P (0) -f gluebaU . (24) 

The reggeon approach cannot tell us anything on the relation between jfa = rn 2 lueball ,t 2 = 
m liuebau an d 7(*i = 0,t 2 = 0). The only thing that we can claim is that the signature factor 
takes into account the steepest part of t-behavior. Therefore, in the next subsection we will 
assume that 

7(*1 = m2 glueballM = m 2 glueball) = 7(*1 = 0, t 2 = 0) | (25) 

this is an extreme assumption which can be used to obtain an upper bound on the cross section. 
Uncertainties related to this and other assumptions we make will be discussed in detail in section 
3.4 and in the summary, section 5. 

3.3 The magnitude of the cross section 

Using Eq. (p2|) ,Eq. (|24]) and Eq. ( p5|) we can rewrite Eq. fl2~0|) in the simple form 



da . ( i i \ 2 G F ^ 2 f s \ 

^U=o = 2vr [a P rn glueball ) x ^ x R l—s j . (26) 

For M H = 100 GeV, the factor S/M 2 H ■ s is equal to 400 GeV 2 for s = 1 GeV 2 . Therefore, 
we can take R e i ~ 0.175 ( see Fig.5 ) for the Tevatron energies. Eq. (|26|) leads to 

j 

— \ yH=Q (M h = 100 GeV, = 1800 GeV) = 6.4pbarn , (27) 
dyu 

This is a very large number, especially if we recall that the total inclusive cross section for 
Higgs production in perturbation theory is on the order of 1 pbarn ]16| . However this estimate 



does not yet contain the suppression due to the (small) probability of the rapidity gap survival, 
which will be discussed in section 4, where we present our final results. Since R e \ grows with 
energy, R e i oc s A we expect that the cross section at the LHC energy is approximately 2 times 
larger than the one in Eq. (E7I). 

3.4 Uncertainties of our estimates 

1. Let us start with the value of R e i. We took it from the experimental data, but we 
nevertheless have two uncertainties associated with it. First, Eq. flT9| ) is written for one Pomeron 
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Figure 5: Experimental data for the ratio R e i(s) = a e i(s) / o~tot{s) ■ so = ^GeV 2 and the log is 
taken on base 10. 



exchange while in experimental data at \fs ~ 20 GeV we have about 30% contamination from 
the secondary Reggeons 0. If we try to extract the one Pomeron exchange from the data, it 
reduces the value of cross section for DPHP by 1.7 times. Therefore, the value for the cross 
section can be about 3.8 pbarn rather than Eq. (|27j). The second uncertainty in evaluation of 
R e i is the value of so; even though s = 1 GeV appears in all phenomenological approaches 
|7], 14], we have no theoretical argument for the value of so- However, since the ratio R e \ in 



Fig. 5 is a rather smooth function of energy we do not expect that the uncertainty in the value 
of So can introduce a large error. 

2. We can take into account also the reactions of Eq. (0) and Eq. (|3|). In Eq. (|26]) we would 
then have to substitute 

a DD (s) 

Rel ► Rd = Rel H , (28) 

O~tot 

where a DD is the cross section of the double diffraction dissociation. Unfortunately, we do not 
have conclusive data on this cross section. However, recent CDF measurements |T7) show that 
this cross section could be rather large ( about 4.7 mb at the Tevatron energy). 

3. The principle uncertainty, however, is associated with the continuation from t = rn 2 glueball 
to t = 0. This is a question which at present can only be addressed in the framework of different 
models. For example, in Veneziano model instead of r/ + ( y t) (see Eq. fll5|) ) a new factor 
appears, namely 

<(t i ) = r(2-a P (t i ))e l ^ M . (29) 

Eq. ( ]2"9"| ) does not give the factor of 7r/2 in Eq. (|24J) and, therefore, decreases the value of the 
cross section given by Eq. ( p6|) by a factor of 2.5. We will return to the discussion of the analytic 
continuation in the summary section. 
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4. As we have discussed in section 2, we can evaluate the value of j(ti = rrig lueball ,t2 = 
m2 giuebaii) = Igiuebaii only if M H /M T < 1. The accuracy of Eq. (||) is 0{Mjj/M^) and we 
thus believe that Eq. Q2"2"D gives a reasonable estimate of igiuebaii for Higgs meson with M H < 
100 GeV. 



4 Survival of large rapidity gaps 

As has been discussed intensively during the past decade (see Refs. |19|, |2l], |22], [23|, |4|, ^5 



2B] , |2"7| , [2"5| , the cross section of Eq. fl2"BD has to be multiplied by a factor S^ pect , which is the 



survival probability of large rapidity gap (LRG) processes. The "experimental" cross section is 
therefore given by 

da(pp -» ppH) 2 do P {pp -» ppH) 

j y ly=o - ^ spect Ty \ y =o ■ m 

Here, dap(pp —> ppH) /dy denotes the cross section calculated in Eq. fl26|). The factor S] pect 
has a very simple meaning - it is a probability of the absence of inelastic interactions of the 
spectators which could produce hadrons inside the LRG. We have rather poor theoretical control 
of the value of the survival probability; this fact reflects the lack of knowledge of the "soft" 
physics stemming from non-perturbative QCD. Different models exist, leading to the values 
about Sg pect ~ lO^ 1 at the Tevatron energies. For double Pomeron processes, this quantity 
has been discussed in Ref |[30| . The result of this analysis is that the value of the survival 
probability for double Pomeron production is almost the same as for "hard" dijet production 



with LRG between them. Fortunately, the value of Si. t has been measured |3l| , and is equal 
to 0.07 for the highest Tevatron energy. 

Multiplying Eq. (p7|) by S 2 spect = 0.07 and taking into account suppression due to the factor 
of Eq. (p9|) , we obtain 

do 

' yH= o(M h = 100 GeV, = 1800 GeV) = 0.2pbarn. (31) 



dy H 



This estimate is not our final result yet, since we still have to correct it by the additional 
suppression factor S par which describes the probability of the absence of the parasite gluon 
emission around the Higgs production vertex (see Fig.4-b) 0. As was argued in Ref.||, 



with 



S par = e -< N G(&y=MM 2 H /s ))> ^ ^2) 

< N G (Ay = HKM) > = ^y~(^ = N^)) „ 2 . 4 . (33) 

N h adrons{one mimjet) 



It gives S^ ar = 0.14 + 0.014. 



The appearance of this factor can be illustrated by the following argument: one of the most 
important differences between the diagrams of Fig. 2 and in Fig. 4 is the fact that the Pomeron 
exchange is almost purely imaginary while the glueball exchange leads to the real amplitude. 
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Imaginary amplitude describes the production of particles and the Pomeron is associated with 
the inelastic process with large multiplicity. Therefore, normally, in a large rapidity region 
Ay = ln(M^/s ) we expect to see a large number of produced particles while in Fig.2 we 
require that only one Higgs boson is produced. Therefore, it seems reasonable to expect a 
suppression for the double-diffractive Higgs production, and this suppression can be described 
by Eq. fl3|) and Eq. @. 

Finally, for the Tevatron energy we expect 

da{pp -> ppH) //-_ 1STl A q2 x9 2 dapjpp^ppH) 
\y=Q\V s — l-oieVJ - o spect x b par x — \ y=0 

= (0.0038 -T- 0.028) pbarn (34) 

Extrapolating to the LHC energy, we have two effects that work in different directions: the 
rise of the Pomeron contribution and the decrease of the Si t with energy. From Ref. [ 30 1 we 
expect that Sg pect (y/s = 16TeV)/S 2 pect (^/s = 1.8 TeV) ~ 0.7 while the rise of the Pomeron 
exchange leads to an extra factor of 2 in Eq. (|34|). Therefore, our final estimate for the LHC is 

dajpp -> ppH) / x 2 2 do P {pp -> ppH) 
^ \y=o(Vs = 16 TeV J = S spect x S 1 ^ x — \ y=0 

= (0.0015 -T- 0.042) pbarn . (35) 



Eq. (|34]) and Eq. (|35]) give significantly larger (by about 5 times) larger cross sections than ex- 
pected for double-diffractive production in pQCD [[33). However, Ref. || contains an estimate 
of the upper bound on double Pomeron Higgs production in pQCD obtained by choosing the 



largest possible value for S par ( see Eq. ([33]) ). This upper bound appears to be about 7 times 



larger than the highest value in Eq. 



5 Summary and discussion 

The approach suggested in this paper is based entirely on non-perturbative QCD. We believe 
that such an approach is logically justified for diffractive Higgs production since even pQCD 
calculations show that this is, to large extent, a "soft" process (see Eq. ([5]) and the following 
discussion). However, just because of this, we have to stress again that the accuracy of our 
calculation is not very good. We feel, however, that our results support the idea that in 
pQCD approach to diffractive Higgs production the running QCD coupling has to be taken at 
the "soft" scale Q 2 ~ 1 GeV 2 . As was argued in Ref. [El, in BLM prescription |52| of taking into 



account the running QCD coupling one can insert the quark bubbles only in the t-channel gluon 
lines in Fig. 3. Therefore, the running QCD coupling depends on the transverse momenta of 
these gluons, and they are determined by the "soft" scaleQ. The Eq. (p~3|) indeed does not depend 
on the QCD coupling, demonstrating the non-perturbative, "soft" character of the discussed 
process. 



1 In this soft regime, the dependence on the coupling constant in the Pomeron can disappear as a consequence 
of scale anomaly p4|. 
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Figure 6: Mueller diagram UJB] for "soft" inclusive Higgs production. 



We obtain quite large values for the cross section of the diffractive Higgs production - after 
integration over the Higgs rapidity y in Eq. ([34|) and Eq. (^) we get 



a(pp -> ppH) (y/s = 1.8 TeV) = 0.019^0.14 pbarn. (36) 

and 

a(pp -> ppH) (V« = 16 TeV") = 0.01^0.27 pbarn. (37) 



Comparing our estimates with the ones based on pQCD [|l], 0, 0, 0, @| we conclude that the 
lowest of our values of the cross section of double Pomeron Higgs production is about the same 
as the highest one in the pQCD approach. However, both our approach and the pQCD one are 
suffering from large uncertainties, stemming from the analytical continuation in our approach 
and from the survival probability of rapidity gap and the absence of "parasite emission" S^ ar 
in pQCD. 

Let us point out that Eq. ( p6| ) shows that the double Pomeron Higgs production constitutes 
a substantial part of the total inclusive Higgs production. Moreover, our calculations lead to 
an additional contribution to the inclusive cross section which is shown in Fig. 6. (Note that 
the triple Pomeron interaction gives a very small contribution to the process in Fig. 6 due to 
the small real part in the Pomeron exchange flSGfl.) Using the same approach as in derivation 



of Eq. fllTD we obtain 

da md (pp ^ H + X) _ n ^ 2 gi (0)g 2 (0)(G p RR ) 2 1 (s 

| VH =o - 7fllti-U,*2-UJ _____ _ | __ 



dy - n(16nB R y A R \M, 



(3* 
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where m 0.5. As a first approximation we can take ( see Eq. (B 

7r i i 2a'r,m 2 f 

7^i = 0,t 2 = 0) = -2^—^ (39) 

where my is the mass of the / - meson which is the first resonance on the secondary Reggeon 
trajectory, and oi R m 2 = 1.5. Substituting Eq. ( j3~9|) in Eq. fl38|) we obtain 

da md (pp^H + X) (G P RR ) 2 B R i 1 / S \ 

s - x ^ x 22Gf s x • (40) 

Eq. (|0D gives 

d* mcl (pp H + X) ^ = 3 410 - 7m6 (41) 

which does not yet contain the suppression arising from the analytical continuation. We take 
this suppression into account by multiplying Eq. (fn|) by factor S 2 ar = 0.14 — 0.014. Unfortu- 
nately, we do not know the value for the ratio G RR /g. In the triple Pomeron parameterization 
of the cross section of diffractive dissociation in hadron reactions [|37j] this ratio changes from 1 
to 0. For G RR /g = 1 we get for the "soft" inclusive cross section the value of 43 -v- 430pbarn. 
On the other hand, taking Field and Fox value |37| for this ratio we obtain a much smaller, but 



still very sizeable value of 0.43 -r- 4.28pbarn. It is thus clear that the evaluation of the "soft" 
contribution to the inclusive Higgs production is plagued by large uncertainties; however, it 
might be bigger than the pQCD one [|1(| . 

We hope that this paper will help to look at diffractive Higgs production from a different 
viewpoint, and will stimulate a much needed further work. To our surprise, despite the very 
different non-perturbative method used here, our estimates for the double-diffract ive produc- 
tion turn out to be not that far from the pQCD calculation |33| ( the average is about 5 times 
larger ). It adds some confidence in both approaches and gives us a hope that one will be able 
to perform a reliable calculation in the nearest future. 
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